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ABSTRACT 

A diversity of resonance configurations may be formed under different migration of two 
giant planets. And the researchers show that the HD 128311 and HD 73526 planetary 
systems are involved in a 2:1 mean motion resonance but not in apsidal corotation, 
because one of the resonance argument circulates over the dynamical evolution. In 
this paper, we investigate potential mechanisms to form the 2:1 librating-circulating 
resonance configuration. 

In the late stage of planetary formation, scattering or colliding among planetesi- 
mals and planetary embryos can frequently occur. Hence, in our model, we consider 
a planetary configuration of two giants together with few terrestrial planets. We find 
that both colliding or scattering events at very early stage of dynamical evolution can 
influence the configurations trapped into resonance. A planet-planet scattering of a 
moderate terrestrial planet, or multiple scattering of smaller planets in a crowded plan- 
etary system can change the resonant configuration. In addition, collision or merging 
can alter the masses and location of the giant planets, which also play an important 
role in shaping the resonant configuration during the dynamical evolution. In this 
sense, the librating-circulating resonance configuration is more likely to form by a 
hybrid mechanism of scattering and collision. 
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1 INTRODUCTION 

To date, more than 40 multiple planetary system have been 
detected beyond our solar system. The commensurability of 
orbital periods is very ubiquitous in the extrasolar plan- 
etary systems. At present, four resonant pair of planets 
(GJ 876, HD 82943, HD 128311, HD 73526) are reported 



to be trapped in 2:1 mea n motion resonanc e (iMarcv et al 



200ll : iMavor et aP l2004l : IVoet et all l2005l : iTinnev et al 
2006h . In recent years, numerous researchers have exten- 



sively investigated the dynami cs and origin of the 2:1 reso- 
nanc e i n the planeta r y syst e ms (iGozdziewski fc Macieiewskil 



2001'; 'Lee & Peale^ '2002'; 'Hadiidemetriou' '200_ 
Hadiidemetriou fc PsY chovo s 200 3: Ji et al. ^ 2003j| b 



nadiiaemetriou ic rsvcnovo s ZUU 6; Ji et al. ^ / 
2OO4I: iKlev. Peitz.fc BrvdenH2004l: iBeauee et al 



2003 



Gavon fc Boij I2OO8I : iLee fc Thomm3 l2009l : IVovatzis et all 
20091 ). The two resonance variables for 2:1 resonance, di = 
Al — 2X2 + TUi , 62 = \\ — 2A2 -I- VJ2 (where A, w are the mean 
longitude and the longitude of periapse respectively, the sub- 
scripts 1, 2 denote the inner and outer planets), are catego- 



rized to librate about: (1)0° in symmetric configuration, (2) 
both 0° and 180° respectively, in the so-called antisymmetric 
configuration, and (3) other d egrees different from 0° or 180° 
in asy mm etric configuration jHadiidemetriou |2002| ; |ji et al.l 
l2003allbll3 : iBeauge et ai]|2003l . l2006l : lLeell2004l ). The GJ 876 
was r evealed to be the first 2:1 resonant system (jMarcv et al.l 
I2OOII ) near an M dwarf star. The GJ 876 system is in apsi- 
dal corotation where the mean motion resonance variables, 
9i and 62 librate about 0° with quite slight amplitudes. On 
the origin of mean motion resonances in the system, a for- 
mation scenario is that they were assembled by migration of 
planets. In the formation of giant planets, if two planets are 
massive enough to open gaps and not far away from each 
other in the disk, the material between the region of them 
can be rapidly cleared off. And then, the dissipation of the 
stuff outside the outer planet and inside the inner planet 
may still force two planets approach each other. Any pro- 
cess that makes two bodies approach each other, which are 
originally separated ap propriately, could result in mean mo- 
tion orbital resona nce (|Lee fc Pealdl2002l: [Kiev et al.ll2005l : 
iMasset et al]l2006l ). 
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Lee and collaborators explored the origin and diversity 
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of the 2:1 mean motion resonance (iLee fc Pealell2002l : lL^ 
I2OO4I : iKlev et all l2005l : iLee et al.l I2OO6I ). Thevset the two 
planets on coplanar, circular orbits with, 02/01 ~ 2/1. The 
outer planet can migrate to ward th e center star at a rea- 
sonable velocity (|Wardll 19971 : see also lLee fc Pealell2002l and 
references therein), 



Table 1. Dynamical orbital fits for two systems. 
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where v = aH^Q, is the kinematic viscosity expressed using 
the Shakura-Sunyaev a-prescription, while H, Q = 2tt/P 
are the height scale of the protoplanetary disk and the or- 
bital angular velocity at a (where P « 27ra^/^(G'Af*)~-^''^), 
and M* is the mass of the center star. In (1), the values of 
Q = 4 X 10""^ and H/a = 0.05 are typical in classical models 
of protoplanetary disks. If we adopt = IM©, the migra- 
tion velocity of the planet with a = lAU, is 9.4 x 10~^ yr~^ 
by (1). They find that such migration of the planets can re- 
sult in resonance capture with the eccentricities of the plan- 
ets growing quickly. If the eccentricity damping induced by 
planet-disk interaction is considered, the system may remain 
stable over secular timescale after the resonance capture has 
happened. In such circumstances, the eccentricities reach 
nearly constant values. The investigations imply that the 2:1 
symmetric resonant configurations may be easily established 
under the forced migration due to planet-disk interaction. 

The numerical explorations and theoretical analysis 
show that the 2:1 resonance planetary configurations could 
be quite div erse due t o different migrat ion in a slightly ec- 
centric disk (|Leell2004l : iKlev et al.ll2005l ). For example, it is 
shown that the orbital solutions of the resonant pairs de- 
rived from Keple r ian fit, could result in unstable behaviors 
l|Vogt et al.ll2005l : iTinnev et al]|2006l ') over the timescale of 
several thousand years in the HD 128311 and HD 73526 
systems. The best-fit dynamical orbital solutions are given 
in Table 1, which each system can remain stable over 10* 
yr. Indeed, both of the systems are stabilized by the 2:1 
resonance, however, with a librating-circulating resonance 
configuration. Planetary configurations such as listed in Ta- 
ble 1 are not observed to lead to the convergent migration 
scen ario mentioned by L ee et al. (2002, 2004, 2006). How- 
ever, iTinnev et al.l (|2006l ) mentioned several possibilities of 
the origin of librating-circulating resonance configuration. 
The resonance configurations may be formed either through 
rapid migration or migration with initial p lanetary eccen- 
triciti es or via a dynamical scattering event, [skndor fc Kiev! 
l|2006l ) provide a mixed evolutionary scenario for such res- 
onance configuration by combining an adiabatic migration 
progress and a sudden gravitational perturbation. In their 
study, firstly, they show the planets could be captured into 
mean motion resonances via inward migration. And then, 
they describe a scenario of a sudden stop of inward migration 
after the resonance capture happened, which is supported 
by o bservati ons of you ng protoplanetary disk (ICalvet et al.l 
l2005i : D'Alcss io et all [2005: Massct ct al. 200^). Finally, a 
small body (~ IOA/q) is scattered during the follow-up 
dynamical evolution of the system. Therefore, the reso- 
nance configuration of the system is turned into a librating- 
circulating stage, and the eccentricities of the planet oscillate 
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with large amplitudes. The results of the numerical studies 
are consistent with the dynamical behaviors expressed using 
solutions listed in Table 1 of HP 128 3 11 and HD 73526 p lan- 
etary systems (|Sandor fc Klevll2006l : fSandor et al.ll2007l '). 

The terrestrial planets are rocky planets, with masses 
ranging in 1 — IQM^. The observations show that 
there are several systems containing several super Earths 
in close-in orbits, e.g., Gl 581 and HD 40307 (see 
|http://exoplanet.eu/|), and recently the nearby solar-like 
star 61 Virginis was reported to harbor two Neptune-like 
planets at about 0.22 AU and 0.48 AU B and a super 
Earth (|Vogt et al.l l2009l). This indicates that terrestrial or 
Neptune-like planets are very common to survive in the 
late stage of planet formation. On the other hand, in 
the numerical studies on the terrestrial planet formation, 
2 — 4 te rrestrial planets are formed with moderate eccen- 
tricities llChambersll200ll : iRavmond eral]|2004 12005| . I2OO6I : 
IZhang fc Jil 20091 ) . It is worthy to pay attention to direct 
or indirect influence upon the resonance configuration in 
the very beginning of evolution. At that time, giant planets 
have been already formed and stopped migrating, while the 
terrestrial planets may have been created but unstable. If 
the periods of the giants are commensurable, the terrestrial 
planets with rather moderate eccentricities could take effect 
on the dyna mical evolution of a system. For exampl e, the 
earlier works (|Sandor fc Klevl2006l : ISandor et al]|2007l ') show 
that librating-circulating configurations of the 2:1 mean mo- 
tion resonance can be formed by combining migration pro- 
cesses and sudden perturbations. In this work, we investi- 
gate potential mechanisms of making a librating-circulating 
resonance configuration. And a novel hybrid mechanism is 
proposed to explain this issue. In the model, we consider the 
planetary configuration of the two giants together with few 



^ The terrestrial planets with masses as small as 2 Earth mass 
and a period of several months orbiting a K dwarf may be de- 
tected in the Habitable zones, by observing 60 groups of 10 nights 
with HARPS on the ESO 3.6 m telescope (M. Mayor 2009, pri- 
vate communication). Ground-based microlensing is sensitive to 
Earth-mass planets orbiting at 1 AU (D. Bennett 2009, private 
communication) around late-M stars, while MPF is sensitive to 
planets down to 0.1 Earth mass in the Habitable zone (near 1 
AU) of solar-type, G and K stars. The future E-VLTs and space 
missions (e.g., GAIA and SIM) will be also hopeful to discover 
multiple Earth-like planets residing in the terrestrial region in a 
planetary system. 
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terrestrial planets. We show that a scattering of a moderate 
mass planet or multiple continuous scattering of the ter- 
restrial bodies can form librating-circulating configuration. 
Moreover, we also find that the collision and merging can 
play a vital role in tuning to a hbrating-circulating mode of 
the giants. In Section 2, we reproduce a scenario of migration 
and capture into a 2:1 resonance. In Section 3, we introduce 
the model and initial setup for our study. We present the 
main results in Section 4. We conclude the outcomes in Sec- 
tion 5. 



2 MIGRATION AND CAPTURED INTO 2:1 
RESONANCE 

As mentioned previously, the 2:1 resonant configurations 
are easily achieved under convergent migration of two gi- 
ant planets. Follow ing N-body damping method given by 
iLee fc Peal3 (120021 '). we modified the MERCURY package 
ijChamberd Il999h l;o examine the situation of migration 
when the star mass is close to IMq. As shown in Figure 
1, the two planets are placed on coplanar circular orbit sim- 
ilar to GJ 876 planetary system, with a2 = 2ai = lAU. 
The migration velocity of the outer planet can be obtained 
by computing the coefficient (~ 5 x 10~^ yr~^) of the 
right side of (1) by substituting Af* for IMq. The ec- 
centricity damping of the outer planet is assumed to be 
62/62 = —100 1 02/02!. The calculations are carried out us- 
ing a hybrid symplectic integrator in MERCURY package. 
When the outer planet is forced to migrate toward the center 
star with 02/02 = —5 x 10~^ yr~^ , the eccentricities of the 
the two planets are excited rapidly, and then they are indeed 
captured into 2:1 mean motion resonance with two resonant 
variables 9i, 62 librating about 0°, while the eccentricities 
equilibrate gradually. The migration scenario is plausible. 
The anti-symmetric, symmetric, and asymmetric configu- 
rations of 2:1 mean motion resonances have been formed, 
which has also been applied to model th e formation of the 
GJ 876 and HD 82943 planetary systems ("Lee & Peale"2002'; 
|Lei|2004; Lcc ct al. 2006; Bcauge_et al 2006) by differen- 
tial migration of planets with constant masse s and initially 
coplan ar nearly-circular orbits. In addition, iBeauge et al.l 
l|2006r ) further pointed out that if they suffer an adiabatic 
process of a slow migration for their orbits, the t wo g ia nt can 
be finally captured into in apsidal corotations. iLed (120041 ') 
showed that some 2:1 resonant configurations may origin 
from an alternation of the planetary mass ratio of two plan- 
ets under migration or multiple-planet scattering in packed 
systems. In the following, we will show that both collision 
and scattering events among the terrestrial bodies and the 
giant planets can yield the 2:1 librating-circulating configu- 
rations. 




3 MODEL AND NUMERICAL SETUP 

In the late stage of planet formation, even though the re- 
gion of the inner disk contain little material after the giants 
formed, several terrestrial planets may still form through 
accreting planetesi mals or embryos in a wider region of 
their feeding zones (IChambersll2(30ll : iRavmond erallbOO^ . 
I2OO5I . I2OO6I : IZhang fc Jill2009l '). Another alternative scenario 
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Figure 1. Time evolution of the semi- major axes, eccentricities, 
and resonant angles of two giant planets in the simulation. The 
initial configuration for the giants was set on coplanar circular or- 
bits reference to GJ 876 system with a2 = lAU and a\ = 0.5AU. 
The migration velocity |a2/a2| is adopted to be 5 X 10~^ yr~^, 
and the eccentricity damping of the outer planet is assumed to be 
62/62 = —100 |a2/a2i. The eccentricities can be excited quickly, 
and then equilibrate gradually after two giants are deeply locked 
into a symmetric 2:1 mean motion resonance, where both of the 
resonant variables d\, 62 librate about 0° with quite small ampli- 
tudes. 



is that the small bodies may be captured into mean mo- 
tion resonances with the giant planets under migration 
(|Malhotra et al.ll2000^ . The shepherd mechanism of the gi- 
ants can also drive a swarm of planetesimals toward the ter- 
restrial region as the sources for further formation. In this 
sense, it is reasonable to assume the system is gas-free, and 
the migration of the giants has been stopped when the ter- 
restrial planets formed. It is supposed that there are a pair 
of giants with commensurable orbits and several terrestrial 
planets inside the region of the inner giant in the system. 
The planetary configuration is much closer to physical sce- 
nario and convenient to study the diversity of resonance sce- 
nario, through dynamical evolution of the terrestrial planets. 

In this work, we cho ose HD 128311 sys tem as the basic 
framework for the model. IVogt et al.l l|2005l ) mentioned that 
the dynamical stable solutions of the HD 128311 are not 
unique. Here we modify the eccentricities of the two giants 
as our initial eccentricities in order to simulate a little ear- 
lier state of the system (see Table 2) . The other parameter s 
for the giant planets are the same as ISandor fc KlevI l|2006l ). 
which are listed in Table 1. The modified orbital parame- 
ters are dynamically stable at least over 10 ^ yr under mu- 
tual g ravitational interaction. In the work of lSandor fc KlevI 
l|2006l ). the semi-major axis of an additional small planet is 
set to be 0.5AU, about half of the semi-major axis of the 
inner planet (~ 1.09AU). As Lee et al. pointed out that a 
small additiona l planet is easy to be captured into mean mo- 
tion resonance (iLee fc PealeH2002l : lLeell2"oo3 ). it indeed hap- 
pened before scattere d for the scenario of the inner planet 
llSandor fc Klevll2006D . In our model, we follow this strategy 
and place 2 — 4 terrestrial planets about 0.5AU to examine 
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Table 2. Initial eccentricities of two giants. 



Ml Inner 


Ml Outer 


M2 Inner 


M2 Outer 


Eccentricity 0.15 


0.46 


0.05 


0.01 



the evolutior0. We consider the condition that there is no 
terrestrial planet or planetary embryo between the orbits 
of two giants in the beginning, indicating that they are en- 
tirely cleared off. The semi-major axes of terrestrial planets 
are distributed in the range 0.3 — 0.75AU. The maximum 
value of the semi-major axis is a little larger than that of 
the 2:1 resonance regime of the inner giant at ~ 0.7AU. The 
terrestrial planets are not well-separated in initial trajec- 
tories, therefore, they can approach each other and further 
merge into a larger body soon after the orbits evolve. 

The eccentricities of the terrestrial planets are ran- 
domly adopted in the range of 0.1 and 0. 3, which is com- 
parable to those of planet a ry formation JChamberd l200ll : 
[Raymond et al.ll2004 l2006l : IZhang fc Jill2009l ). Certain ini- 
tial eccentricities mean that the planetary orbits may in- 
tersect in the beginning of evolution. The planetary masses 
also play an important role in the colliding or scattering 
event. Consequently, the terrestrial planets occupy a total 
mass of 10 — ISM®, and these bodies are placed inside the 
inner giant. As regards to the specific mass of an individual 
terrestrial planet, we treat the issue by following the route 
of unequal-mass or equal-mass. In the former case, the mass 
of a terrestrial body is varied randomly among 3 — 6M0, 
and the total mass is required to meet the above condition. 
For equal-mass, the planetary mass is SM^ , and , 

respectively, corresponding to the number of 2, 3 and 4 for 
the terrestrial planets in the system. In this instance, the 
overall mass is lOM^ or 12M(g. Moreover, the density of the 
giants and terrestrial planet is, respectively, considered with 
respect to that of Jupiter and Earth in solar system. Hence, 
we can obtain the size of each body, and the collision and 
merging can be also taken into account in the run. In each 
simulation, the terrestrial planets are in coplanar orbits ref- 
erence to the giant planets. The other angles of the mean 
anomaly M and the longitudes of periapse cc are randomly 
distributed between 0° and 360° for each orbit. 

We prepare 120 simulation s for two group s and use 
the hybrid symplectic integrator (|Chamberslll999l ) in MER- 
CURY package to integrate all the orbits. The algorithm is 
capable of handling close encounters among the bodies. In 
the calculations, we consider that close encounters will take 
place while the bodies are separated by not more than their 
3 Hill radii. And the collision and merging occurs when the 

As aforementioned, iRavmond et al.l | |2004|) show that a few ter- 
restrial planets with masses close to the Earth's mass can be 
formed within 1 AU. However, due to the limitation of current 
observational precision, the Earth-like planets similar to those 
in solar system residing in the terrestrial region have not been 
discovered yet. It is likely that all the super Earths recently dis- 
covered may be formed farther from the central stars and then 
migrated to their present positions. This may imply the difficulty 
in making two or more of such planets in nearby orbits remain 
stable during the evolution. The adopted hypothesis should be ex- 
amined by future high-precision observations, to unveil multiple 
stable terrestrial planets in the Habitable zones. 



minimum distance between any of the two objects is equal to 
or less than the summation of their physical radii. In addi- 
tion, we adopt 4 days as length of time step, which is about 
one twentieth period of the terrestrial body residing at 0.3 
AU. Most of the simulations are carried out for 10'^ — 10^ 
yr. Nevertheless, a few simulations are performed to extend 
on longer integral time. 



4 RESULTS 

The previous works (|Sandor fc KlevI I2OO6I : ISandor et al.l 
I2OO7I I on librating-circulating configurations of the 2:1 mean 
motion resonance in the HD 128311 and HD 73526 systems 
have reported that such resonance configuration can be pro- 
duced by scattering of an additional small planet through 
orbital migration of a sudden stop. Our results also confirm 
their conclusion. Moreover, we further show that collision 
and merging can play a major role in dynamical evolu- 
tion and both mechanisms can modify the resonance an- 
gles of the giants. The configuration engaged in the mean 
motion resonances is varied or destructed. Concretely, in 
our results, there are 16 simulations involved in inducing 
librating-circulating resonance configuration, among which 
about 56% are shaped by mixed events of planet-planet 
scattering and merging (or collision). In this section, we 
present major results of several simulations that illustrate 
the dynamical interaction of planet-planet scattering and 
collisions. In Table 3 are listed the initial conditions of ter- 
restrial planets in four runs. Next we will address the possi- 
bility to form librating-circulating configuration with certain 
initials for terrestrial planets during dynamical evolution. 

4.1 Planet-planet scattering 

The present observations show that the eccentricities of the 
extrasolar planets can be much higher on average than those 
of the solar system. Many potential mechanisms (such as 
secular perturbations from a massive planetary companion 
or passing stars, interactions of orbital migration via mean 
motion resonance, or planet-plan et scattering in the plan- 
etary systems, etc., see details in lFord fc Rasidl2008l ) have 
been proposed to explain the origin of the distribution of 
high eccentric i ties. B ase d on the realistic initial conditions, 
iFord fc Rasiol (|2008l ') and lChatteriee et al.l (|2008l ) show that 
the combination of planet-planet scattering and tidal circu- 
larization may lead to form several close-in giants and re- 
produce the observed eccentricity distribution. In addition, 
in the following we show that planet-planet scattering can 
also shape the librating-circulating configuration trapped in 
resonance. 

Librating-circulating configuration of mean motion res- 
onance may be form ed by a single plane t scattering event 
with moderate mass (|Sandor fc Klevll2006l '). Herein, we show 

^ In the following, we refer to collision as a merging event be- 
tween a small terrestrial planet and a giant planet, while merging 
specifically, means the circumstance between two small terrestrial 
planets. MERCURY package models two planets inelastically to 
merge together and produce a single new body by conserving mass 
and total momentum. 
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Figure 2. (a) Time evolution of the mass of T2, and tlie semi-major axes of Tl, T2, T3 for simulation M2120. (b)Time evolution of the 
resonance angle 0i, semi- major axis and eccentricity of T2 in simulation M2120. T3 merges into T2 at about one hundred yr, because of 
their very close initial orbits. After Tl merges into T2, the eccentricity of 'Merger' T2 is greatly excited by gravitational perturbation of 
the inner giant, and then T2 is scattered at ~ 3 X 10* yr. The configuration is modulated to be a librating-circulating stage after such 
scattering event. 



Table 3. Initial conditions of the terrestrial planets for 4 runs. 
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M{deg) 


M1204 Tl 


5.0 


0.66 


0.14 


344 


226 


T2 


5.0 


0.63 


0.29 


171 


238 


M1205 Tl 


5.0 


0.54 


0.15 


92 


219 


T2 


5.0 


0.39 


0.11 


280 


225 


M1227 Tl 


3.0 


0.48 


0.12 


51 


132 


T2 


3.0 


0.34 


0.17 


63 


203 


T3 


3.0 


0.54 


0.27 


56 


209 


T4 


3.0 


0.58 


0.16 


118 


245 


M2120 Tl 


3.7 


0.39 


0.12 
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39 


T2 


5.9 


0.57 


0.10 


144 


32 


T3 


4.9 


0.56 


0.15 


344 
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" We use four digits following the letter M to label each simu- 
lation. The first digit denotes the models using different initial 
eccentricities of the giant planets, as shown in Table 2. The next 
digit indicates an alternative strategy as treating the individual 
masses of the terrestrial planets in the simulation (the label 1, 2 
respectively for unequal and equal masses). The last two digits 
are used to stand for each simulation. The letter T means terres- 
trial planets. In each run, we set all terrestrial planets on coplanar 
orbits with respect to those of the giant planets. 

an interesting example that is consistent with such mecha- 
nism but related to a more complicated physical scenario. 
Figure 2(a) illustrates the time evolution of the mass of T2, 
and the semi-major axes of three additional terrestrial plan- 
ets in simulation M2120. From Table 3, we can clearly notice 
that T2 and T3 are set close in the initial setup, and then 
they are rapidly merged together into a larger body at about 
100 yr. Due to significant excitation of the eccentricities, the 
orbital crossings of the terrestrial planets quickly occur and 
subsequently cause the merging event between T2 and T3. 
In Figure 2(b), we can see the 'Merger' T2 remains in a 
stable orbit for about 10* yr but its eccentricity oscillates 
with larger amplitudes. And then the consequence of stir- 



ring the eccentricity is that T2 is scattered at ~ 3 x 10* yr. 
Here comes the librating-circulating configuration of the pair 
giants, where the resonance argument 6\ modulates with 
moderate amplitudes about 0°, and 62 circulates as time 
evolves. The difference is that, in our assumptions, the ini- 
tial conditions of the two giants are simply required to meet 
aijax ~ 2/1 (where the subscripts 1, 2 denote the inner and 
outer giant, respectively, hereinafter), so their configuration 
after migration is not required to be exactly trapped into 2:1 
resonance, where two resonant angles are initially both in 
circulating mode. In our model, we are not concerned about 
how two giants planets move to the orbits close to the 2:1 
resonance through certain consequences of migration. But 
the initial setup we require is that several terrestrial planets 
stay inside the region of the inner giant in the gas-free disk 
after the termination of migration for two giants. 

We find that if the masses of the scattered terrestrial 
planets are low, then the changes of the resonance configura- 
tions may reinstate after some while. Generally, small planet 
scattering is capable of increasing the amplitude of the res- 
onance angles. It is very likely that continuous scattering 
events resulting from small planets can alter the configura- 
tion of the mean motion resonance. Figure 3 illustrates the 
results of the simulation of M1204. In this run, there are two 
terrestrial planets with equal masses of mri = mT2 = SM® 
initially moving about the giants. From Table 3, we know 
that, though their semi-major axes are originally approxi- 
mate to each other, the two terrestrial planets can escape 
colliding from each other. Because the orbit of Tl is a bit 
broader and less eccentric than that of T2, and the difference 
in the mean longitudes of periapse is AA = Ai — A2 = 161° 
(where \ = w -\- M, M is the mean anomaly), near 180°. 
As Figure 3 shows, the 2:1 resonance of two giants is ini- 
tially in symmetric configuration where {61,62) ~ (0°,Q°). 
At the time the first terrestrial planet Tl is scattered at 366 
yr (marked by triangles), the amplitudes of the resonance 
angles become enlarged, but the apsidal corotation is main- 
tained. Until the second terrestrial planet is scattered at 593 
yr, soon after the first scattering event, the apsidal corota- 
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Figure 3. Evolution of tiie 2:1 mean motion resonance variables 
6*1, 02, and the eccentricities of the two giants ei, 62, for the sim- 
ulation M1204. The sign triangles label the time of two terrestrial 
objects that scattering happens. Obviously, the resonance config- 
uration after two sequent scattering events of terrestrial planets 
turns into a librating-circulating state from an original librating- 
librating mode of {6*1, 6I2) « (0°,0°). 

tion is broken up. And the eccentricities of two giants ei and 
62 fluctuate with large amplitudes, which is similar to the 
HP f28311 and H P 73526 systems (jSandor fc KlevI l2006l : 
ISandor et ai]|2007l V In the circumstances, our outcomes are 
well consistent with those of ISandor fc Kiev! (|2006l '). 

4.2 Merging and Collision 

As shown previously, the librating-circulating configuration 
can arise from a planet-planet scattering of a moderate mass 
planet or multiple continuous scattering for small terrestrial 
planets in the various initial conditions. In this section, we 
will further focus on the mixture role of merging and collision 
to generate a librating-circulating mode in the dynamical 
evolution. 

Figure 4 shows the results of M1227, which are involved 
in the very complicated process of colliding and scattering 
events. In the initial stage, there are 4 terrestrial planets in 
total, and each terrestrial planet bears a pretty smaller mass 
of 3M0. In this scenario, it is more likely to observe mul- 
tiple scattering or colliding events (marked by upper letter 
'S', 'C, respectively in Fig.4) over timescale of the dynam- 
ical evolution. From the figure, we can see that the inner 
giant is hit by T3, and then T4 is scattered within 500 yr. 
However, we may notice that the 2:1 mean motion resonance 
configuration are not significantly changed, where the first 
planet-planet scattering event or collision with the inner gi- 
ant before 500 yr do not modify the oscillating status of two 
resonance arguments 9i and 62, due to their smaller plane- 
tary masses {3M^). In the subsequent evolution, Tl and T2 
are merged together to form a larger terrestrial body, and the 
eccentricity of the merger is gradually excited to high values, 
which triggers the collision between a new 6M^ body and 
the outer giant. The second collision event of the merged ter- 
restrial planet and the outer giant occurs at about 3500 yr. 
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Figure 4. Dynamical evolution of the 2:1 mean motion resonance 
variables 8i and 62, and eccentricity eyi semi-major axis a^i of 
the terrestrial planet Tl, and the eccentricities ei, 62 of the two 
giants for the run of M1227. The upper letters 'C and 'S' denote 
a Collision or Scattering event occurs at the corresponding time. 
Collision and scattering events are observed in the simulation. 
However, the first scattering event or collision with the inner giant 
before 500 yr cannot change the vibration status of the resonance 
angles. The configuration of mean motion resonance has been 
triggered to be a librating-circulating mode at 3461 yr, because 
Tl and T2 are merged into one new bigger terrestrial body, which 
collides with the outer giant planet at about 3500 yr. 



which induce the two resonant angles 9i and O2 to go from 
librating-librating phase into librating-circulating phase. It 
is worthy to note that the eccentricities of the giants un- 
dergo smaller oscillations. The reason for the variations of 
the eccentricities is the conservation of the total momen- 
tum, while the impactor loses the momentum and transfers 
to the outer giant, which gains the momentum. Then the 
eccentricity 62 is slightly modulated if the semi-major axis 
remains and further influences ei of the inner giant. 

Figure 5 illustrates the evolution of simulation M1205. 
It is similar to the case of M1227. However, the difference 
is that the first 'C denotes a merging event between Tl 
and T2 at about 500 yr. It is not difficult to understand 
that such merging scenario for two terrestrial planets can- 
not alter the masses of giant planets, and will not cause 
the variations of their semi-major axes and eccentricities. 
However, subsequent evolution is to be expected soon af- 
ter the first collision. The second collision occurs between 
T12 (as a merged larger body) and the inner giant at about 
1300 yr. In this case, two resonant arguments 9i and 82 vary 
from a librating-librating phase into the librating-circulating 
state. Similarly, both of the eccentricities for two giants are 
slightly modulated. The above two figures show that a ter- 
restrial planet's colliding with either an inner or outer giant 
may finally shape librating-circulating configuration of mean 
motion resonances during the dynamical evolution. 

In addition, there are lots of simulations that show the 
destruction of the systems, where one of the two giants is 
scattered faraway to interstellar space or even entirely out 
of a planetary system. In most of our runs, we notice that 
both scattering and colliding events may coexist in the evo- 
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Figure 5. Time Evolution of tiie simulation M1205. The first 
upper letter 'C denotes a merging event between Tl and T2. 
And the second collision occurs between T12 (as a merged body) 
and the inner giant. It is clear that colliding with inner giant can 
also form a librating-circulating configuration engaged in mean 
motion resonances. 

lution. The planet-pla net scattering model l|Ford fc Rasid 
I2OO8I : IChatteriee et al.l 2008) can well explain the distribu- 
tion of relatively higher eccentricities of the extrasolar plan- 
ets in comparison with the observations, and further indicate 
that much more planets can be excited in the eccentricity 
in the dynamical evolution. In the late stage of planet for- 
mation, the evolution of orbits is quite chaotic, because at 
that time the terrestrial planets could frequently intersect 
with each orbit. In this sense, it is very likely for terrestrial 
planets or massive planetary embryos to suffer from gravita- 
tional scattering and collisions in the system. In final, such 
hybrid mechanism may participate in determining and form- 
ing a 2:1 librating-circulating resonant configuration in the 
planetary system. 

In order to examine the limiting conditions for the 
above mechanism to work, we further performed over 100 
sets of simulations to find the minimum masses of the 
terrestrial planets, which convert librating-circulating sta- 
tus for the giant planets. As for planet-planet scattering 
case, in simulations related to M1204, we find that when 
mri = mT2 ~ 3.8M(g, the first terrestrial planet Tl hit 
the inner giant at about 128 yr, and the second terrestrial 
planet is scattered at 269 yr soon after the first collision 
event, and the librating-circulating resonance configuration 
is induced. In other simulations of this group, the terres- 
trial planet masses mri = mT2 < S.SMq cannot trigger the 
librating-circulating resonance configuration. As for merg- 
ing and collision scenario, the additional runs of M1205 
shows that when rriTi = mr2 = 2.5M0, a merging event 
between Tl and T2 occurred at about 368 yr, then the 
second collision occurs between T12 (as a merged larger 
body) and the inner giant at about 671 yr, finally result- 
ing in the librating-circulating configuration of mean motion 
resonance. Also the simulations related to M1205 show that 
the librating-circulating resonance configuration of two giant 
planets cannot be formed, when the terrestrial planet masses 
rriTi = mT2 < 2.5Af®. It is not difficult to understand that 
the presented minimum masses of terrestrial planets have 
different values in these simulations, to alter the status of 



the mean motion resonance configuration of the giants. This 
is because such conversion of the resonant planetary configu- 
ration may mainly not only depend on the terrestrial planet 
masses, but also rely on their other initial conditions as well. 



5 CONCLUSIONS AND DISCUSSIONS 

In this work, we have investigated potential mechanisms to 
shape 2:1 librating-circulating resonance configuration by 
considering a planetary system of two giants accompanying 
with few terrestrial planets with coplanar orbits in the late 
stage of planet formation. In the model, the system is con- 
sidered to be gas-free, and two giants with commensurable 
orbits stop migrating and several terrestrial planets have 
formed inside the inner giant. The configuration is much 
closer to real planet ary formation scenario according to nu- 
mcrical simulations (|Chambersll200ll : iRavmond et al]|2004l : 
I Zhang fc J i 20091). In conclusion, we summarize the main 
results as follows. 

In the late stage of planetary formation, planet-planet 
scattering or colliding among planetesimals and embryos can 
frequently occur. Our results show that not only a single 
planetary scattering of a terrestrial planet with a moder- 
ate mass can result in the librating-circulating configura- 
tion, but several continuous planetary scattering with rather 
smaller terrestrial masses is also at work. Additionally, if 
two giant planets are initially engaged in a 2:1 symmetric 
resonance and their eccentricities oscillate with large am- 
plitudes, the collisions arising from the giants and other 
small bodies may change librating amplitudes of the reso- 
nance angles during the evolution. If the apsidal corotation 
is disintegrated, the configuration may turn into a librating- 
circulating status. Obviously, the more mass of a perturbing 
terrestrial body may have much greater infiuence on the 
commensurable giant planets. In most simulations, colliding 
and scattering events can be found and they can increase or 
decrease the fluctuation in the amplitude of the resonant an- 
gles, even dramatically destroy the whole system. In a word, 
the librating-circulating configuration of mean motion reso- 
nance is likely to generate by a mixed mechanism of colliding 
and scattering. 

In addition, the librating-circulating configuration 
could be generated through long-term evolution of planetary 
formation. A librating-circul ating config uration trapped in 
3:2 mean motion resonance (|Zhang &: Jil |2009| ) is unveiled 
during the formation of the terrestrial planets, which is sim- 
ilar to the newly-d iscovered planetary system HD 45346 
(jCorreia et al.l [20091 ). They show that the two terrestrial 
planets were formed within 50 Myr (see their Figure 6 for 
details), and after that time the frequent orbital crossings 
of them and their interaction with the inner giant planet 
may finally lead to a capture of a 3:2 resonance. Such reso- 
nant configuration is believed to hold over hundreds of Myrs. 
However, more resonant configurations for less massive plan- 
ets are expected to reveal by future higher accuracy space- 
based projects in search for terrestrial planets (such as TPF, 
Darwin, SIM). The innovative findings will encourage one to 
more carefully study their dynamics and origin. 
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